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ABSTRACT 



> 

' Context. The detection of radio pulses from cosmic ray air showers is a potentially powerful new detection mechanism for studying spectrum and 
composition of ultra high energy cosmic rays that needs to be understood in greater detail. The radiation consists in large part of geosynchrotron 
radiation. The intensity of this radiation depends, among other factors, on the energy of the primary particle and the angle of the shower axis with 
^~>) respect to the geomagnetic field. 

J — Aims. Since the radiation mechanism is based on particle acceleration, the atmospheric electric field can play an important role. Especially inside 
^— ^ thunderclouds large electric fields can be present. In this paper we examine the contribution of an electric field to the emission mechanism 
■»««^ theoretically and experimentally. 

Methods. Two mechanisms of amplification of radio emission are considered: the acceleration radiation of the shower particles and the radiation 
Oh from the current that is produced by ionization electrons moving in the electric field. For both mechanisms analytical estimates are made of their 
effects on the radio pulse height. We selected lopes data recorded during thunderstorms, periods of heavy cloudiness and periods of cloudless 
weather. We tested whether the correlations with geomagnetic angle and primary energy vary with atmospheric conditions. 
Results. We find that during thunderstorms the radio emission can be strongly enhanced. The present data suggests that the observed amplification 
^ is caused by acceleration of the shower electrons and positrons. In the near future, extensions of lopes and the construction of lofar will help to 
. . identify the mechanism in more detail. No amplified pulses were found during periods of cloudless sky or heavy cloudiness, suggesting that the 
^ electric field effect for radio air shower measurements can be safely ignored during non-thunderstorm conditions. 
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1. Introduction in the atmosphere. Furthermore. IWilsmil (Il957h and ICharmanI 

(1 19681) suggested that ionization electrons, left behind by the 

The first detection of radio pulses from extensive air show- g^S, emit radiation when accelerated in a background electric 

ers (HAS) was in 1964 d Jelley et alj [19651) and several emis- ^gj^ ^^^^^^ thunderclouds the electric field can be large enough 

sionmechanisms have been proposed to explain them.lAskaryanI accelerate the ionization electrons up to energies high enough 

(11962!) calculated the Cherenkov radiation resulting from the produce ongoing ionization (typically 100 kVm '). This effect 

negative charge excess in extensive air showers.| Kahn & Lerche| ^^Ugd runaway brea kdown and a ca lculation of die associated 

(11966) considered two more mechanisms, both driven by the ge- ^^^^^ p^j^g ^j^^g bv lGurevich et alJ (H^. 
omagnetic field. Firstly, the geomagnetic field separates the neg- 
ative and positive charges. The electric dipole created in this way 

will emit Cherenkov-like radiation in the atmosphere. Secondly, !„ the 1970's several groups did experiments with shower 

the transverse cuiTent that is generated by the charge separation arrays and radio antennas. One of their aims was to infer the 

produces a radiation field which is strongly beamed forwards in correct emission mechanism from polarization measurements, 

the direction of the EAS. They predicted the latter mechanism to Although most experiments were in favour of the transverse 

be dominant. current mechanism, a large contribution of other mechanisms 

The atmosp heric elec t ric fie ld also contributes to the total could not be excluded. Large spreads in radio intensity, inabil- 

radio emission. ICharmanI d 19671) calculated the effect of charge ity to filter out radio interference, and the significant but un- 

separation due to electric fields and concluded that this effect known effects of atmospheric conditions led to the abandonment 

can contribute significantiy when a large electric field is present of these experiment (see Baggio et al. 1977). Excessively large 

radio pulses of EAS during thunde rstorms were found experi- 

Send offprint requests to: S.Buitink@astro.ru.nl mentally bv iMandolesi et al.l (|1974|) . 
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Recently, iGurevich et alj (|2004|) reported the discovery of 
radio pulses of duration ~ 0.5 fis associated with EAS during 
thunderstorms at the Tien Shan Scientific Station in Kazakhstan. 
These were detected with antennas sensitive to frequencies be- 
tween 0. 1 and 30 MHz. The forte satellite has detected strong 
intracloud radio pulses in the 26-48 MHz band, which are 
correlated with discharge processes inside the thundercloud 
jjacobson 2003J. 

The develo pment of lofar rev i ved th e interest in radio de- 
tection of EAS. iFalcke & GorhamI (|2003|) describe the emission 
in terms of coherent geosynchrotron emission. Although no ex- 
plicit comparisons are made, they expect this mechanism to be 
largely equivalent to the transverse current mechanism, since 
it also results from charge separation and is beamed forward. 
Detailed simulati ons of the geosynchrotron emissi on from EAS 
are p resented in iHuege &~F alcke (2003) and Hue ge & Falckd 
(l2005h . 

The LOPES experiment was set up at the kascade array site 
in Karlsruhe, consisting of 10 and later 30 radio antennas sen- 
sitive to radiation in the band 40-80 MHz. The presence of the 
KASCADE particle detectors on the site, new techniques and mod- 
ern hardware allow higher resol ution and a highe r detection rate 
than the old experiments. In falcke et al.l (l2005 l) it was shown 
that the measured antenna electric field is strongly correlated 
with the muon number and the angle of the shower axis with 
the geomagnetic field. The first correlation is a strong indication 
that the radio emission from EAS is coherent, while the latter 
correlation proves that the dominant emission mechanism is of 
geomagnetic nature. 

In this paper we investigate theoretically and experimentally 
the effect of atmospheric conditions on the radio emission. We 
compare sets of events that were recorded by lopes under vari- 
ous weather conditions: fair weather, large nimbostratus clouds, 
and thunderstorms. Under violent weather conditions, the effect 
of geoelectric mechanisms increases. When this increase is large 
enough it may dominate over the geomagnetic emission mecha- 
nism. Understanding the eff'ect of atmospheric conditions on the 
emitted radiation is crucial for a correct determination of the pri- 
mary energy from observed radio pulses. 

2. Electric fields inside clouds 

We present some general characteristics of thunderstorms based 
on the information provided by MacGorman & Rust ( 1998). In 
fair weather, i.e. atmospheric conditions in which electrified 
clouds are absent, there is a downward electric field present 
with a field strength of ~ 100 Vm ' at ground level. The field 
strength decreases rapidly with altitude and has values below 
10 Vm ' at altitudes of a few hunderd meter and higher. The 
associated fair weather current charges up cloud boundaries, be- 
cause clouds have lower conductivity than the free atmosphere. 
Other effects, such as ion capture and collisions between po- 
larized cloud particles, also contribute to the charging up of 
clouds. Clouds can typically gain field strengths of a few hun- 
dred Vm"'. Nimbostratus clouds, which have a typical thickness 
of more than 2000 m can have fields of the order of lOkVm 
The largest electric fields are found inside thunderstorms, where 
locally field strengths can reach values up to 100 kVm"' . In most 
clouds this field is directed vertically (either upwards or down- 
wards, depending on the type of cloud), but thunderclouds con- 
tain complex charge distributions and can have local fields in 
any direction. Thunderclouds can have a vertical extent of ~ 10 
km. The electric field at ground level is strongly affected by the 
electric processes inside thunderclouds. Although it can not be 



used to estimate the field strength inside the cloud, a change in 
the (polarity of) the ground field is a strong indication that large 
electric fields are present overhead. In the context of EAS radio 
emission, ground level electric field mills can be used as a warn- 
ing system for violent electric phenomena in the atmosphere. 



3. Electric field influence on radiation 

The atmospheric electric field acts on the radio emission from 
EAS in various ways. We distinguish two generations of elec- 
trons: the relativistic electrons from pair creation in the EAS 
(called shower electrons from here) and the non-relativistic elec- 
trons resulting from the ionization of air molecules by the EAS 
particles (called ionization electrons from here). The shower 
electrons are created together with an equal amount of positrons. 
In this section, we will use Gaussian units. 

1. The electric force accelerates the shower electrons and 
positrons, producing radiation in more or less the same way 
as the magnetic field does. This effect is described in Sec. 
EH 

2. The ionization electrons are accelerated in the electric field. 
A radio pulse will be emitted from the short current that is 
produced in this way. This eff'ect is described in Sec. 13.21 

3. As the electrons and positrons move through the electric field 
they can gain or lose energy, which has an influence on the 
electromagnetic cascade. When the electric field points up- 
wards and the shower is vertical the electrons gain as much 
energy as the positrons lose. However, when a charge excess 
of e exists a shower with primary energy £0 = 10'* eV and 
electron number ~ 10^ moving through a thundercloud 
sustaining a 1 kV/cm vertical electric field over 2 km, the 
energy gain of the whole shower will be: 



AE = eNq\E\Ax ^ 2 ■ 10"*eV 



(1) 



where l^l is the electric field strength, and we have used a 
typical value of e = 0.08. This corresponds to ~2% of the 
primary energy. This value is only a very rough approxima- 
tion. 

4. The electric field acts on all charged particles in the EAS and 
can thus influence the longitudinal and lateral development 
of the EAS. This may aff'ect the coherence and the shape and 
size of the radio footprint. 

The last two effects can best be investigated by includi ng elec- 
tric field effects in a Monte Carlo code like corsika (Hec k et all 
1998) and are not studied in this paper. The first two effects will 
be further explored below. Although both mechanisms can be re- 
sponsible for a strong enhancement of radio emission, they can 
be distinguished by their temporal and spectral radiation profiles. 

3.1. Acceleration of air shower electrons 

The radiation part of the electric fiel d of a moving electric charge 
can be expressed as (lJacksonlll975h : 



E{x, t) = 



(2) 



where e is the unit charge, R is the distance to the observer, 
P - v/c is the velocity of the charge and « is a unit vector point- 
ing in the direction of the observer (approximating the index of 
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refraction for radio waves in the atmosphere to be unity). The 
associated vector potential is: 



A(x, t) 



(^) 



1/2 



[RE], 



(3) 



As the electron and positron paths are curved in the geomag- 
netic field they emit synchrotron radiation as in Fig. [1] From 
Eqns. |2]and|3]the vector potential in frequency domain can be 
calculated (Huege & Falck&.2003l) : 



A(R,co) = 



a>e 



.'(w£-f) I 



-e||A||(w) ± ej.Ai(w)] 



(4) 



where the plus-sign corresponds to electrons and the minus-sign 
to positrons, <y is the angular frequency of the radiation, the 
unit vector in the plane of curvature of the trajectories and ej_ the 
unit vector perpendicular to that plane, both transverse to n. The 
components are given by: 



V3c\r 



Ax = e- 



2p 



1/2 



with 



3c 



3/2 



(5) 



(6) 



(7) 



where Ka is the modified Bessel-function of order a, 6 is the 
(small) angle between n and /3, and the radius of curvature p is 
given by: 



P = 



vym^c 
eB sin a 



(8) 



where v and y are the speed and Lorentz factor of the emitting 
particle, B is the magnetic field strength and a is the angle of the 
trajectory with the magnetic field direction (the pitch angle). It 
can be seen from Eqn. |4]that although the radio emission com- 
ponents of the electron and positron perpendicular to the plane 
of curvature cancel out, the components in the plane of curvature 
add up. 



The radiated power of an accelerated charge can be ex- 
pressed in the following form, known as the Lienard result 
(Jackson 1975): 

2 

P^^-y^[{pf-{p>^pf] (9) 

from which it can be seen that the power is related to the acceler- 
ation as P|| oc in case of linear acceleration and oc 7^)8 
when the acceleration is perpendicular to the direction of mo- 
tion. Furthermore, oc y^^Fy and oc y^^F^, where F\\ is a 
force in the direction of motion and Fj^ is a force perpendicu- 
lar to the direction of motion. Comparing the eff'ect of these two 
forces one finds: 

1 „2 



Pn = 



while: 

2 



2 1-2 



(10) 



(11) 



which is a factor y^ greater The Lorentz force is always perpen- 
dicular to the particles' direction, but the electric force can have 
any angle with respect to the trajectory. 




Fig. 2. The electrons and positrons make curved trajectories in 
the magnetic field. Under influence of a downward directed elec- 
tric field, the positrons (electrons) are accelerated (decelerated). 
The asymmetry in the trajectories will be reflected in the radio 



ermssion. 




Fig. 1. The electrons and positrons travel along curved trajecto- 
ries in the magnetic field, emitting synchrotron radiation. 



For a more detailed view we turn our attention to the vec- 
tor potential. We will treat a general case in which a linear force 
and a perpendicular force act on a charge and evaluate the vector 
potentials caused by these two forces. We define the unit vector 
ei as perpendicular to the particle's velocity j8 and lying in the 
orbital plane of the particle. From Eqn.|4]we know that the emis- 
sion from an electron/positron pair is strongly polarized in this 
direction. We therefore evaluate the vector potentials of linear 
and perpendicular acceleration along this unit vector for an ob- 
server in the orbital plane. Since both the charge and in Eqn. 
|2] have different signs for electrons and positrons, the resulting 
field is the same for both types of particles. For linear accelera- 
tion Eqns.|2]and|3]reduce to: 



A0,t) 



jAnc 



sin 6* 



(1 -f5cos6) 



(12) 
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Fig. 3. Electron/positron pairs are created in a horizontal 
shower The downward electric force works in the same direc- 
tion as the Lorentz force for both species. The radio emission is 
amplified. 



and for perpendicular acceleration: 

cos 6»-y8 



(1 -ySCOSfl) 



(13) 



when n lies in the orbital plane of the particle. Note: the sub- 
scripts to A refer to the type of acceleration (linear or perpendic- 
ular) that cause the vector potential. For small 6 these equations 
can be written in the form: 



A||(0,f) = 



r 



yd 



(1 +y202): 



2^ (l+y202)3^^ 



(14) 



(15) 



where the relativistic approximation (1 -yS)"' = 2y^ is used. We 
will now compare the peak values of A\\ and Aj^. The peak of Ay 
lies atyO-l jl and the peak of Aj. at 6* = 0. We find that: 



A||max r'ii\\ ^ y^Fii 



(16) 

(17) 



When an electric field is absent the peak value of the vector po- 
tential is: 



(18) 



When an electric field E with pitch angle 6 (defined as the angle 
between E and j8) is present the peak value of the vector potential 
is: 



^E,max 



oc y^q(cB sin ff ± £■ sin 6) + y^qE cos 6 



(19) 



where the plus sign corresponds to the case where the electric 
force and the Lorentz force act in the same direction and the mi- 
nus sign to the case where they act in opposite directions. The 
Lorentz factor y can diff'er from yo as a result of linear accelera- 
tion or deceleration. Now we can define an amplification factor 

^ amo as . 



An.max r \ L Esmd 

A^amp = = — 1 + ^— 

Ao.max \yo) \ Be sm a 



yo \yo 



E COS 6 



Be sin a 

Three effects can be distinguished that cause amplification: 



(20) 



- The part of the electric field that is directed perpendicular 
to the particles. The amplification depends on the force as 
A^amp ^ F. For a typical case of a horizontal shower in a 
vertical electric field (see Fig. [3]l with values B - 0.5 G, 
a - 25° and yo = 30 the amplification factor is 1.15 for an 
electric field of 1 kVm"' and 17 for a field of 100 kVm"'. 

- The part of the electric field that acts linearly on the parti- 
cles' trajectory. The amplification depends on the force as 
A^amp ^ F/yo- For the case of a vertical shower in a vertical 
electric field and the same characteristic values as mentioned 
above the amplification factor is 1.5 for an electric field of 
100 kVm-i. 

- As the particles are accelerated/decelerated their Lorentz 
factor increases/decreases. The amplification depends on 
the Lorentz factor as A^amp iyljof' ■ This eff'ect depends 
strongly on the track lengths of the particles. Suppose a pair 
is created in a vertical shower in a vertical electric field and 
the same characteristic values mentioned above. When the 
particles have crossed an altitude difference of 100 m the 
positron has reached a Lorentz factor of 50, while the elec- 
tron's Lorentz factor is down to 10 (see Fig. |2]i. The cor- 
responding amplification is (50^' + 10-')/2 ■ 30-^ ~ 2.3. After 
200 m the electron has changed direction and its radiation no 
longer reaches the ground. The positron now has a Lorentz 
factor of 70, corresponding to an amplification factor of 6.3. 
At a height of 4 km the mean free path length of electrons 
and positrons is in the order of a few hundred meter 



Amplification of emission from a complete shower 

So far, only the emission of a single particle pair was discussed. 
Eqn.|20]does not apply to complete showers, consisting of many 
electron-positron pairs with varying pitch angles, energies and 
track lengths. An observer at ground level can see the emis- 
sion of a particle only for a fraction of its lifetime, since the 
particle follows a curved trajectory. When the emission from 
a particle is amplified due to an electric force, this does not 
necessarily mean that an observer will see an increase in emis- 
sion. The observed pulses of single particles have a duration of 
A/ = nplcy^ ~ 10~" s while the total pulse of a shower is of the 
order AT" = L/{2cy^) ~ 10^^ s, where L is the length of the to- 
tal shower. The pulses of individual particles are distributed over 
the period Ar, so the emission is not totally coherent. Since the 
number of particles is much higher than AT /At the emission is 
also not completely incoherent. 

The radius of curvature of a relativistic particle is given by 
Eqn. [8] The time width of a single pulse, as measured by a 
ground observer, depends on the Lorentz factor and the applied 
force as Af oc y'^F]_^. Some of the amplification effects men- 
tioned in the section above wiU therefore not contribute to a total 
increase of emission. 



- The part of the electric force that is directed perpendicular to 
the particle orbit increases the peak value of A proportional 
to Fj_. The time width of the observed pulse goes down as 
Fj} and the contribution to the integrated emission of the 
whole shower remains the same. 

- The part of the electric force that is directed along the parti- 
cle orbit does not influence the time width of the pulse (ini- 
tially). However, the sign of Ay is dependent on the viewing 
angle (see Eqn.[T3]). When the pulses are distributed over AT" 
pulses of opposite polarity will be added and partiaUy cancel 
out. 



Buitink et al: Radio emission of air showers in thunderstorms 



5 



- When the particle is accelerated its field strength increases 
proportional to y^, while the pulse width decreases as y^^. 
The contribution to the integrated shower field strength 
roughly increases with y. The power of the integrated pulse 
increases with y^. 

From this treatment it appears that the increase in particle ve- 
locity is the most important factor in pulse amplification for an 
integrated shower However, when the time width of a pulse de- 
creases because the particles' trajectory is bent over a larger an- 
gle this also means more observers will be able to see emission 
from that particle. In general, an observer will be able to see a 
larger fraction of the total amount of particles and, in effect, an 
amplification of the radio pulse. 

The total amplification of a radio pulse from a complete 
shower depends on the distribution of particles and the position 
of the ground observer and cannot be easily predicted. A de- 
tailed Monte Carlo simulation with realistic particle distribution 
will give more reliable results and is subject of further research. 
For now, we can regard the amplification factor as given in Eqn. 
l20las an upper limit for the observed amplification. 

3.2. Acceleration of ionization electrons 

The shower particles ionize air molecules and leave behind free 
electrons and positive ions. The electrons can recombine with 
the ions in a time scale of seconds, but on a much shorter time 
scale of a few tens of nanoseconds th e electrons atta ch to oxy- 
gen molecules forming negative ions (IWilsoniri957h . When an 
electric field is present the free electrons are accelerated produc- 
ing a current pulse. Because of freq uent c ollisions they gain a 
drift velocity of ~ 100 m/s (Charman 19681) . The duration of the 
pulse depends on the attachment time Tan of electrons to oxygen 
molecules (which is a function of the electron energy) and the 
angle under which the shower is viewed. In contrast to the ra- 
dio emission from the shower electrons/positrons the radiation 
is not beamed forward since the free electrons do not become 
relativistic. The associated frequency up to which the emission 
is coherent is ~ 10 MHz for an observer in the direction of the 
shower (inverse of Tan ~ 100 ns). 

When ionized electrons gain an energy of e > e^. a; 0. 1 - 1 
MeV they can ionize new molecules. If the electric field is strong 
enough to accelerate ionization electr ons to such energi es a pro- 
cess called runaway breakdown (Gur evich et al.|[T992h can oc- 
cur The critical field strength of ~ 100 - 150 kV/m, needed 
for this effect, is present only inside thunderclouds. In the run- 
away breakdown process two generations of electrons are cre- 
ated: relativistic runaway electrons and slow thermal electrons, 
which can both contribute to the radio signal. In Fig.|4]the con- 
tribution of the ionization electrons is schematically displayed. 
A simple current would produce a dipolar radiation field. The 
real radiation pattern will be more complex because of the finite 
length of the current, the transverse width of the current and the 
existence of relativistic electrons. The radiation pattern o f the 
runaway breakdown is calculated in lGurevich et alj (l2002h for a 
vertical shower and resembles that of a current pulse. The pulse 
amplitude is calculated to be several orders of magnitude higher 
than the geosynchrotron emission from the HAS. 

3.3. Distinguishing between the mechanisms 

Both mechanisms can be responsible for an amplification of the 
radio pulse from EAS. There are several ways to distinguish be- 
tween them: 




Fig. 4. An EAS passes through a thunderstorm cloud emitting 
the usual geosynchrotron emission. The ionization electrons are 
accelerated inside the cloud. In this picture the radiation pattern 
is displayed as dipole radiation. The real pattern of a current 
pulse suffers from diffraction. 



1. The radio pulse from the shower electrons will have ap- 
proximately the same width as produced by the geosyn- 
chrotron mechanism and is dependent on the longitudinal 
separation of the air shower particles. The pulse width of 
the ionization current pulse is determined by the electron at- 
tachment timescale (Af^n ~ 100 ns) and the timescale asso- 
ciated with the scale /a o f the runaway breakd own current 
Af ~ IJc * 200 - 300 ns dOurevich et al.ll2004l) . 

2. The pulse from the shower particles will be polarized in the 
plane of curvature, as in the geosynchrotron mechanism. The 
polarization of the ionization current pulse is in the vertical 
plane of the current and the observer 

3. The radiation from the ionization current is emitted in all 
directions, while the radiation from the relativistic shower 
particles is beamed forward, as shown in Fig.|4] 

4. Experimental Setup 

In 2004, the lopes array consisted of 10 dipole ant ennas, placed 
on th e same location as the kascade experiment ( lAntoni et al.l 
|2003|) which provides triggers for lopes and records the muon 
and electron components of the EAS, used to reconstruct the en- 
ergy of the primary and its direction. Presently, lopes measures 
only the polarization in the east-west plane. A layout of the ex- 
periment is shown in Fig.|5] Details about the experiment al setup 
and th e reduction of the data can be found in Horneffe r et al.l 
(2004). Additional information about the weather was obtained 
from a weather station at Kai'lsruhe (49°02'N 8°22'E) in the 
archive of a free weather servefl 

Four sets of data were selected from the 2004 database of 
lopes: 

1. Events with the highest kascade particle count. The events 
in this selection have either a truncated muon number above 
2-10^ or an electron number exceeding 5 ■ 10^ (412 events 
spread over the period January-September). 



' http : //meteo . inf ospace . ru/wcarch/html/index . sht 
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Antenna Layout 



Table 1. Applied cuts on selections. 
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Fig. 5. Layout of the LOPES experiment. Boxes indicate posi- 
tions of KASCADE detector stations, filling up an area of 200 m 
by 200 m, and circles indicate antenna positions 



2. Fair weather events which took place during periods with 0% 
cloud coverage (9455 events spread over the period March- 
September). 

3. Events which took place while the sky was covered by nim- 
bostratus clouds for more than 90% (2659 events spread over 
the period January-March). 

4. Events during thunderstorms, which were identified by look- 
ing at lightning strike map^ and the dynamic spectra of 
LOPES. The radio emission of actual lightning strikes show 
up on these spectra as bright lines because the antenna sig- 
nal is saturated. In our selection we only used thunderstorms 
that were visible both on the maps and the spectra. The radio 
events that were recorded between these strikes or half an 
hour before the first or after the last strike on the spectra are 
regarded as thunderstorm events (3510 events taken from 1 1 
thunderstorms in the period May-August). 

Together, all these events form only a very small fraction 
of the total lopes database, because the weather information is 
not complete and even if it was, most weather conditions do 
not match the criteria of one of these selections. The selections 
include events for which a radio signal was not detected. The 
weather at the lopes site is expected to differ only slightly from 
the weather as archived at the Karlsruhe weather station. 

To determine whether the radio peak is significant or not a 
cross-correlation beam was created. A combined signal of all 
LOPES antennas is reconstructed by temporal shifting of the pulses 
in accordance with the arrival direction in the kascade EAS data. 
The CC-beam is then calculated by adding correlations of all 
possible antenna pairs. The radio signal as a function of time is 
fitted with a Gaussian and is considered a detection when the fit- 
ted peak is larger than the background noise by 3 sigma. Since 
the amplitude calibration of the antennas is not yet completed, 
the pulse heights are given in arbitrary units. In the scope of this 
paper this is not a problem, since we investigate relative differ- 
ences between sets of events. 
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5. Results 



In lFalcke et al.l (|2005[) it was shown that the strength of the radio 
signal depends on the geomagnetic angle as ( 1 - cos a) when 
it is normalized with the (truncated) muon number. We show 
this correlation here, using the selection of events that is the first 
set listed in Sec. |4] and applying the cuts Usted in Table [1] un- 
der 'largest'. In Fig.|6]the normalized pulse heights are plotted 
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Fig. 6. Normalized pulse height of the events from selection 1 
of Sec. |4] plotted against the geomagnetic angle. Excess pulse 
heights are defined as the normalized pulse height minus this fit. 



against geomagnetic angle and a fit is made. The normalization 
is done by dividing by the truncated muon number and multi- 
plying by 10^. In order to compare data points to this fit we will 
calculate the difference between the points and the fit, not the 
ratio. The reason for this is that the ratio is a value that is nor- 
malized for geomagnetic value. For an amplified pulse caused 
by an electric field this normalization is not suitable. The total 
radiation can be seen as consisting of a geomagnetic part (which 
scales with geomagnetic angle) and an electric field part (which 
does not scale with geomagnetic angle). We therefore define the 
pulse height excess as the normalized pulse height minus the fit 
value (the difference between a data point and the fit in Fig.|6]l. 
The error-weighed mean pulse height excess is calculated as: 



1 V 

iyil=0.33 



(21) 



where = 2 o", Y, are the pulse height excess values and cr,- 
the errors. Various effects are included in the calculation of this 



2 2,2,2 
erf = crt;„ + cr% + cr 



http : //webcam . paanstra . nl/ 



bg-^" fit "r " phase <^ge (22) 

where cTbg is the root mean square of the background signal and 
(Tfi, is the error in the Gaussian fit to the measured radio pulse. 



Buitink et al: Radio emission of air showers in thunderstorms 



7 



Table 2. Mean pulse height excess. 



cloudless 


nimbostratus 


thunderstorm 


X 0.17 


-0.15 


1.27 



An error in the phase caHbration of the antennas translates into 
an error a phase in the formed CC-beam. The antenna gain factor 
depends on the signal direction, so an error in the signal direc- 
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Fig. 7. Pulse height excess is plotted against zenith angle for 
conservative cut. 



and thunderstorm events, three cuts were applied, which are all 
listed in Table[T] In the 'conservative' cut the maximum distance 
of the shower core to the array centre is 91 m and the zenith an- 
gle of the shower must be smaller than 50°. Beyond these limits 
the values of the general sho wer parameters may not be well re- 
constructed by KASCADE (see lAntoni et al.l (|2003|) ). Furthermore, 
at least 6 out of 10 antennas must have detected a radio sig- 
nal and the reconstructed muon number must exceed 10^. In 
the 'distant' cut the constraint on distance is dropped. The con- 
servative cut yields only 3 thunderstorms events, 3 fair weather 
events and 2 nimbostratus events. Their pulse height excesses are 
plotted against zenith angle in Fig.|2l One of the three thunder- 
storm events shows a large excess, while the other events deviate 
from the fit within Icr. Because of the low statistics we will fo- 
cus on the results of the 'distant' cut, which leaves 14 out of 
3510 thunderstorm events (0.40%), 15 out of 9455 (0.16%) fair 
weather events and 7 out of 2659 (0.26%) nimbostratus events. 
The detection rate increases considerably during thunderstorm. 
The detection rate during nimbostratus conditions also seems to 
be slightly higher, but one should note that this is all low number 
statistics. 

In Fig. [8] the excess pulse height for events of different 
weather selections are plotted against geomagnetic angle, where 
the distant cut is applied. Figs.l9l [T0l[m and[T2lcontain the same 
data points, now plotted against respectively HAS energy as es- 
timated by KASCADE, zenith angle, azimuth angle and mean dis- 
tance of the antennas to the shower axis. There is a bias towards 
events with a positive pulse height excess, because these have 



a higher chance to be detected on a 3cr level. The mean pulse 
height excesses are listed in Table|2l 



U 



tion translates into an additional error in the gain factor, o-g^.. The ^ 

last error is very small for zenith angles below 50°. For most ra- ^ 

dio pulses (Jhg is dominant, but for strong pulses cr phase gives .5^ 
the largest contribution. All but one data point in Fig. |6]have an 

excess < 2. In the evaluation of the fair weather, nimbostratus ^ 
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Fig. 8. Pulse height excess is plotted against geomagnetic angle 
for distant cut. 
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Fig. 9. Pulse height excess is plotted against shower energy for 
distant cut. 



Some events have extremely nice radio pulses like in Fig. 
[T3] The signal-to-noise ratio in the cross-correlation beam (see 
Fig.[T4li is much larger than for other events and the coherence is 
very high. A selection of these events was made by eye ('radio' 
cut in Table [U and the pulse height excesses are plotted in Fig. 
[TS] All of these events have signal-to-noise ratios of > 10. In 
the fair weather selection, 5 out of 9455 events had such strong 
radio emission (0.05%), in the nimbostratus selection 1 out of 
2659 (0.04%) and in the thunderstorm selection 11 out of 3510 
events (0.3%). The events that show a bright pulse because of the 
shower size appear near the bottom of the plot, because of the 
normalization with muon number Bright pulses from showers 
with a relatively small muon number appear in the upper part of 
the plot. Only thunderstorm events are present in this region. 

To further check the uniqueness of the high excess thunder- 
storm events, new selections of twin events were made from the 
LOPES database. For each thunderstorm event (A through K in 
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Fig. 10. Pulse height excess is plotted against zenith angle for 
distant cut. 
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Fig. 12. Pulse height excess is plotted against mean distance to 
the shower axis for distant cut. kascade reconstruction of muon 
number becomes unreliable above 91m. 
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Fig. 11. Pulse height excess is plotted against azimuth angle to 
the shower axis for distant cut. 
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Fig. 13. Signals of ten lopes antennas for a very radio bright 
event (event K in Fig. \T5[ . (Field strength values are not cah- 
brated.) 



Fig. [T5] l a selection was made of events with approximately the 
same muon and electron number (both within 5%). The zenith 
angle and mean distance to the shower axis of these twin events 
is not necessarily the same. The pulse height excesses, corrected 
for geomagnetic angle by dividing by the fit in Fig. |6] of all 
events are plotted in Fig.[T7]by group. For groups A through F 
the pulse height excess of the thunderstorm event is significantly 
larger than those of their twins, while in groups G through K 
the thunderstorm events have excesses similar to those of their 
twins. The pulse heights are not normalized for mean distance 
to the shower core, since this dependence is not yet clearly es- 
tablished. Any reasonable normalization (e.g. oc exp(-7?/100 
m)) will not change the appearance of Fig. [T7] significantly, i.e. 
the thunderstorm events of groups A through F still have much 
larger pulse height excesses than their twins. 

All the events that show a large excess in Fig.[T5]have signif- 
icantly larger pulse heights than their twins, while the low excess 
events have pulse heights similar to their twins. In group E two 
of the twin events have a large pulse height excess. Lightning 
maps and data from the weather station show that these events 
have also occured under thunderstorm conditions. (They were 



not in the original selection because no lightning strikes were 
visible in the dynamic spectrum during these thunderstorms.) 



6. Discussion 

It is found that during thunderstorms the radio signal from HAS 
can be strongly amplified. Due to the low number of events for 
which the weather conditions could be reliably reconstructed, it 
was necessary to include events that were more than 91m away 
from the array core. The kascade reconstruction of the muon 
number is not fully reliable for showers with large zenith angle 
or a distant core. When the estimation of the muon number is too 
low the calculated excess values will be too high. The conserva- 
tive cut is the most reliable but leaves the fewest data points. 
The distant cut leaves more data points simply because there are 
many events with a distant core, not because they are detected 
more efficiently. Although the number of data points in the con- 
servative cut is too low for statistics, it reflects the structure of the 
plots where the distant cut is applied. Large excesses for thun- 
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Fig. 14. Cross correlation beam of the ten signals of Fig. [T3] 
The blue line is the cross correlated signal and the purple line is 
a Gaussian fit. (Field strength values are not calibrated.) 
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Fig. 16. Pulse height excess is plotted against the geomagnetic 
angle for radio cut. 



Thunderstorm events & Twins 




Zenith Angle (degrees) 

Fig. 15. Pulse height excess is plotted against the zenith angle 
for radio cut. The thunderstorm events are labelled A through K. 
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Fig. 17. Pulse height excess, corrected for geomagnetic angle, 
for thunderstorm events A through K and their twins. In group 
E, two of the twin events are probably also thunderstorm events 
(see text for details). 



derstorm events were found, while the spread for fair weather 
and nimbostratus events is small. 

Although some thunderstorm events show a large excess, 
others fall inside the same spread as fair weather and nimbostra- 
tus events. This can be because the selected time windows for 
the thunderstorm selection probably contains some time before 
and after the thunderstorms. Also, even when the thunderstorm 
is at its strongest, the amplification of the radio emission will de- 
pend on the local field distribution inside the cloud and the angle 
of the shower axis with the electric field. 

During nimbostratus conditions no amplified radio emission 
was found. This could be because the electric field is too weak, 
or because these clouds have a smaller vertical extent than thun- 
derstorm clouds and most showers reach their maximum above 
the nimbostratus cloud, where no large electric field is present. 

Figs. |9] through [T2l show no correlations between the pulse 
height excess and other EAS parameters (primary energy, zenith 
angle, azimuth angle and distance to the shower axis), which in- 
dicates that the observed amplification is caused by the weather 
condition. Fig. [16] and perhaps Fig. [8] seem to indicate that the 
amplified pulses cluster at low geomagnetic angles. The statis- 



tics are not good enough to make a statement on whether this 
feature is coincidental or not. More data is needed to study this 
possible correlation. 

The twin events in Fig.[T7]have the same muon and electron 
number, but other EAS parameters vary. The values in this plot 
are normalized with muon number and corrected for geomag- 
netic angle. The spread in pulse height excess within a group of 
twin events can be due to differences in zenith angle, azimuth 
angle or distance. The spread is, in groups A through F, small 
compared to the pulse height excess value of the thunderstorm 
event(s). When, in the future, a larger database of events is avail- 
able a similar analysis can be done for twin events that also share 
the same zenith angle, azimuth angle and distance. 

Radio pulses observed with lopes typically have a width of 
50 - 60 ns. Due to the 40-80 MHz band filtering any radio pulse 
shorter than that will appear broadened. Broader radio pulses, 
however, do maintain more or less their original width. The ob- 
served pulses of the amplified thunderstorm events have widths 
of ~ 50 ns (see e.g. Fig.O like all other lopes events, and are 
probably not ionization current pulses, which would have widths 
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of at least 100 ns (see Sec. 13.31 ) We suspect therefore that it is the 
direct influence of the electric field on the shower electrons and 
positrons that is responsible for the amplification of the radio 
emission during thunderstorm conditions. 

The small distances between the antennas does not allow an 
evaluation of the lateral distribution of the radio emission and 
LOPES only measures one polarization, so no additional tests can 
be done to identify the mechanism at the moment. 

7. Conclusion 

It is shown that during thunderstorm conditions the radio emis- 
sion of EAS is largely amplified. There are two mechanisms 
which can explain this amplification: acceleration radiation from 
the shower electrons and radiation from a current pulse of (run- 
away) ionization electrons. The measured pulse widths (~ 50 ns) 
suggest that the latter cannot be the observed mechanism. To 
identify the mechanism with more certainty, more information 
about the real pulse width, lateral distribution and polarization is 
needed. At the moment, the lopes experiment is unable to pro- 
vide this information, but with future additions to the experi- 
ment, such as dual polarization, it will be possible to resolve this 
problem. Also, in the short future, lofar stations will be able 
to help find the answer, since they occupy a larger ground area, 
operate in a wider frequency range and measure polarization. 

For both the lopes and the lofar experiment it is advisable to 
keep detailed weather information, like cloud coverage, ground 
level electric field and the occurence of lightning strikes. 

With LOFAR it will be possible to trace lightning activity by 
three dimensional imaging. This technique allows localization 
and mapping of lightning strikes, and possibly also thunder- 
storm processes emitting weaker radiation such as stepped lead- 
ers and high altitude lightning. This off'ers unique opportunities 
and promises significant further advances in this area. 
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